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p^ t Abstract. We analyse the extension of Chiral Perturbation Theory to describe a me- 

son gas out of thermal equilibrium. For that purpose, we let the pion decay constant be 
a time-dependent function and work within the Schwinger-Keldysh contour technique. 
A useful connection with curved space-time QFT allows to consistently renormalise the 
model, introducing two new low-energy constants in the chiral limit. We discuss the 

iy~\ ■ applicability of our approach within a Relativistic Heavy-Ion Collision environment. 

In particular, we investigate the formation of Disoriented Chiral Condensate domains 
in this model, via the parametric resonance mechanism. 

ON 

ON 

^ INTRODUCTION AND MOTIVATION 

The forthcoming experiments on Relativistic Heavy Ion Collision (RHIC) at BNL 
and CERN will be able to test accurately the dynamics of the QCD plasma. After 
the collision, the plasma formed in the central rapidity region cools down via hydro- 
dynamic expansion, and nonequilibrium effects become important in that regime. 
Among them, one of the most interesting suggestions has been the formation of 
the so called Disoriented Chiral Condensates (DCC). The DCC were proposed 
originally in [1] as misaligned vacuum regions, where the chiral field points out in 
directions in isospin space different from that where the vacuum expectation value 
of the pion field vanishes. If such regions were formed, one could observe large 
clusters of pions emitted coherently from the plasma as the pion field relaxes to 
the normal vacuum. This kind of behaviour is indeed observed in Centauro and 
anti-Centauro events in cosmic ray experiments [2]. However, one should point out 
that a clear signal for DCC formation has not been observed yet in the RHIC at 
Fermilab [3], although it could well happen that the DCC's are too small to be 
directly detected and one has to think of alternative observables (see below). 
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On the other hand, after the hadronisation time, a proper description of the mi- 
croscopic meson dynamics makes it compulsory to use an effective low-energy theory 
for QCD. In this context, the chiral symmetry plays a fundamental role. The effec- 
tive theory must incorporate all the QCD symmetries and the chiral spontaneous 
symmetry breaking (SSB) pattern, so that the Nambu-Goldstone bosons (NGB) 
are the lightest mesons (it, K, rj). The light quark masses are then introduced 
perturbatively. One possible choice is simply an O(N) model, with the standard 
classical SSB potential. Its fundamental fields are the N — 1 pions and the a, 
which acquires a nonzero vacuum expectation value v. This is the Linear Sigma 
Model (LSM) description. However, one should bear in mind that the LSM be- 
comes nonperturbative in the coupling constant at low energies, so that alternative 
perturbative expansions have to be used, such as large N. Besides, the LSM only 
shares the QCD chiral symmetry breaking pattern for N — 4. An alternative ap- 
proach is to construct an effective theory as an infinite sum of terms with increasing 
number of derivatives, only for the NGB fields. This is the description based in the 
Nonlinear Sigma Model (NLSM), which is the lowest order action one can write 
down in this expansion. Higher order corrections come both from NGB loops and 
higher order lagrangians and can be renormalised order by order in energies, yield- 
ing finite predictions for the meson observables. The unknown coefficients, which 
encode all the information on the underlying theory, absorb the loop infinities and 
can be fitted to experiment. This framework constitutes the so called chiral pertur- 
bation theory (ChPT) [4,5]. The perturbative expansion is carried out in terms of 
the ratio of the 0(p) meson energy scales of the theory (masses, external momenta, 
temperature and so on) and the chiral scale A x ~ 1 GeV (see [6-9] for a review). 

Nonequilibrium effects such as the DCC's have been investigated in the literature 
using O(N) models with initial thermal equilibrium conditions cr(t — 0) = and 
7T a (t = 0) = 0. In this context, two different scenarios for DCC formation have 
been proposed: the first one takes place in the early stages of the plasma evolution. 
Roughly speaking, as the field rolls down along the potential, long wavelength 
modes grow exponentially and enhance the formation of DCC's. There have been 
several approaches in the literature to implement this idea in the O(N) model, like 
classical simulations [10], large N [11], or analysis based on reasonable assumptions 
on the kinematics [12,13]. Typical DCC sizes within this approach are of the order 
of 2-3 fm, containing n n ~ 0.2 fm 3 pions, whereas the plasma cools down in a 
proper time of about r ~ 5-10 fm/c. As commented above, these numbers yield 
too small DCC's to be observed directly. A second suggestion, which has been 
recently proposed [14-16] is based on the parametric resonance mechanism and 
inherits the idea from inflationary reheating [17]. In this approach, the a field is 
very close to the bottom of the potential but it is still oscillating around it (it 
clearly overshoots the vacuum if the initial conditions are imposed on the top) in 
the late stage of the plasma evolution. Those oscillations transfer energy to the 
pion modes, giving rise to exponentially growing pion solutions for certain bands in 
momentum space, via parametric resonance. Recent work within the LSM in this 
approach predicts rather larger DCC, of sizes up to 5 fm [15]. More details about 



this mechanism will be given below. 

In the present work, we will explore the applicability of ChPT to describe the 
meson plasma out of thermal equilibrium. So far, this formalism has been applied 
only in equilibrium, to study the low T (T = 0(p)) meson gas and the chiral phase 
transition [18,19]. The key idea of our approach is to make use of the derivative 
expansion consistently defined in ChPT in order to study the system not far from 
equilibrium. It is therefore best suited for the late stage evolution and has the 
additional advantages typical of standard ChPT, i.e, it deals only with NGB fields 
and is equally applicable to three flavours. We will show that a systematic power 
counting can be established in this case and, furthermore, that the renormalisation 
program can be consistently implemented. Details can be found in [20] . In addition, 
in the last section we will explore the possibility of describing DCC formation within 
ChPT, via parametric resonance. 

THE MODEL AND CHIRAL POWER COUNTING 

Our starting point is the nonlinear sigma model (NLSM) where we let the pion 
decay constant- the only relevant parameter to the lowest order in derivatives- be 
time dependent. In the context of a RHIC, such time dependence can be thought 
of as proper time evolution within the so called Bjorken initial conditions [21], 
where observables depend only on proper time and not on rapidity. This picture 
is consistent with the experimental observations. We take the initial time t = 0, 
having in mind that it would correspond to a proper time r ~ 1-2 fm/c, a typical 
hadronisation time. Thus, we will consider the following NLSM action 

S[U] = f dt f d 3 x f —p- tr d^(x, t)d»U(x, t) (1) 

Here, C is the Schwinger-Keldysh contour (see [20] for details), which 
parametrises the nonequilibrium path integral where we are considering thermal 
equilibrium for t < at a temperature Tj = /3" 1 , as the initial condition. Note 
that the action (1) is chiral invariant (U — > LUFtf) by construction, which will play 
an important role in what follows. As a first approximation, we will be interested 
only in the strict chiral limit for two light flavours, i.e, massless pions. Therefore, 
we are not including any explicit symmetry-breaking term in the action. Thus, 
f(t < 0) = / ~ 93 MeV to leading order (/ ^ /„. to higher orders) and for t > 
the system departs from equilibrium. Note that, since we choose that departure to 
be instantaneous, f(t) cannot be analytical at t — 0. This is just an artifact of the 
approximation and should not have any effect on the long-time behaviour. Finally, 
as customary, U(x) is parametrised in terms of pion fields n a as: 

u(x, t) = jL | [f {t) _ ^ V2 j + irX j (2) 

and 7i a (ti — ifii) = 7i a (ti) is the equilibrium boundary condition, with tj < 0. 



The new ingredient we need to incorporate in the power counting in order to be 
consistent with ChPT is then 
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and so on. Obviously, our results will depend upon the choice of f(t). One can 
think of /(£) as an external source, to which we wish to obtain the nonequilibrium 
response of the system. Alternatively, this model can be thought of to lowest 
order as the LSM with the time-dependent constraint a 2 + n 2 = f 2 {t). We shall 
discuss below a reasonable assumption for f(t) in connection with DCC formation. 
Meanwhile, we shall keep /(£) arbitrary. 

To lowest order in the pion fields, the above NLSM action can be written as 

5 [tt] = -\ [ d 4 x7i a (x, t) \u + m 2 (t)l iT a (x , t) (4) 



2Jc 

where J c d 4 x = J c dtJ d 3 x and m 2 {t) = —f(t)/f(t). That is, the model accom- 
modates a time-dependent pion mass term, without breaking explicitly the chiral 
symmetry. This effect is the same as switching on an external curved space-time 
background, as we will see in the next section. 

RENORMALISATION AND CURVED SPACE-TIME 

Once we have defined our nonequilibrium power counting, we can apply ChPT 
to calculate the time evolution of the observables. In doing so, we must pay special 
attention to renormalisation. The fact that there is a time-dependent mass term 
indicates that there can be new time-dependent infinities in the chiral loops. How- 
ever, we are in the chiral limit, so we are not allowed to introduce the usual 0{p 4 ) 
mass and wave function counterterms breaking the chiral symmetry [5]. In other 
words, we should be able to construct the most general fourth order action, which 
in particular should include new terms (and hence new low-energy constants) to 
cancel those extra divergences, preserving exactly the chiral symmetry. 

In order to find this 0(p A ) lagrangian, we will make use of a very fruitful analogy: 
the action (1) is equivalent to formulate the NLSM on a curved space-time back- 
ground corresponding to a spatially flat Robertson- Walker metric, with scale factor 
a(t) = f(t)/f(0 + ) (see [20] for details). Note that in this language, m 2 {t) in (4) 
represents the minimal coupling with the RW metric preserving chiral invariance. 

Therefore, we can construct the 0(p A ) action as: 

S 4 [U, g,R]= f d A x^Tg [C A [U, g] + (L n Rg^ + L l2 R^)tid^d v U] (5) 

where g is the metric determinant, C±\U,g\ stands for the standard (equilibrium) 
lagrangian [5] with indices raised and lowered with the g^ v metric and the rest are 



new 0(p A ) invariant couplings with the scalar curvature R(x) and the Ricci tensor 
Rfiu{x) in the chiral limit. These are the new terms we need, where L\\ and L\ 2 
are the new coupling constants. In fact, this problem has been already considered 
in [22] in order to study the energy-momentum tensor of QCD at low energies. In 
that work it has been found that L\\ is renormalised in dimensional regularisation, 
whereas L 12 is already finite. Their numerical values can be obtained from the 
experimental information on the QCD energy-momentum form factors. They yield 
L 12 — —2.7 x 10~ 3 and L\ x {p, = lGeV) ~ 1.4 x 10 -3 where \x is the renormalisation 
scale. In our case, with our RW metric we get to 0(7r 2 ), 
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The above lagrangian should take care of the nonequilibrium infinities we might 
find in the pion two-point function. We will see below that this is indeed the case. 



THE PION DECAY FUNCTIONS f n (T) 

The first observable one might think of calculating in ChPT is the pion decay con- 
stant to one loop. In the nonequilibrium model, it will become a time-dependent 
function f n (t). One should point out that the definition of /„- is subtle even in 
thermal equilibrium [19,23]. In addition, one has in general f%(T) 7^ ft{T) cor- 
responding to the axial current spatial and temporal components and due to the 
loss of Lorentz covariance in the thermal bath [24]. We refer to [20] for details on 
how to define properly f n (t) out of equilibrium. Once this has been done, one has 
to consider the one loop diagrams for the pion two-point function coming from (1) 
plus the tree level ones from (6). The final result up to 0(p 4 ) reads [20] 
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for t > 0, with fi^it) in (7) and Go(t) is nothing but the equal-time pion two-point 
function Go{t) = Goix,x) with Goix,y) the solution of the differential equation 



{d x + m 2 (x )} G (x, y) = -5cix° - y°)8®(x - y) 



(10) 



with KMS equilibrium conditions G^(x, i, — i0i',y) = G^(x, t^y), Gq and Gq 
being advanced and retarded correlation functions. Clearly, this equation cannot 
be solved analytically for an arbitrary /(£), but it can be managed numerically. 
Therefore, one must remember that Go(t) depends implicitly on f(t) through (10). 
As a consistency check, the results (8)-(9) reproduce the equilibrium result [25] 
when we switch off the time derivatives of f(t): 

\an 2 = [fi{T)\ 2 = f [\ - ?pj (ii) 

An interesting consequence of our result is that f%(t) ^ fl(t) to one-loop, unlike 
the equilibrium case. In addition, from (8)-(9) and (7) we see that the difference 
[/*(t)] 2 — [flit)} 2 is finite, so that /*(£) and /£(£) can be renormalised at the same 
time, which is another consistency check. We remark that Go(t) contains in general 
UV divergences, to be absorbed by fi(t) and J2{t) in the renormalisation of L n . 
An explicit check of this renormalisation procedure will follow in the next section. 



DISORIENTED CHIRAL CONDENSATES IN CHPT 

In this section we will consider a particular choice of f(t) and apply our previous 
results. Our motivation is the possibility of generating DCC-like structures in this 
context. We shall sketch some of our preliminary results here, while details of the 
calculation and further work will be postponed to a forthcoming paper. 

As we have discussed above, our approach is meant to be useful in a stage of the 
plasma evolution where the departure from equilibrium is of the same order as the 
meson energies. Hence, we should be able to obtain similar results as the analysis 
performed in the LSM in the parametric resonance regime [14-16], where the rolling 
down of the a field is in its late oscillatory period. This is the same behaviour of 
the inflaton field in reheating [17]. One then allows for a time-dependent classical 
background cr(t) in the LSM, splitting the field as 

a(x,t)=a(t) + 5a(x,t) (12) 

where 5a is the quantum fluctuation. As a first approximation, one can neglect 
the pion fluctuations (n 2 ) <C v 2 [14,16] and solve the equation of motion, which 
yields just a(t) = aocosm a t. Here, a is the initial field amplitude, which in this 
approximation is a small quantity. Even though, one can still produce exponentially 
growing pion fields (DCC) which in the end will be responsible for the damping of 
the oscillations as the field relaxes to equilibrium. One should bear in mind that 
neglecting (it 2 ) to lowest order is a rather crude approximation, as pointed out 
in [15], which is clearly not valid for large times when the pion correlator grows 
significantly. Nonetheless, we will carry on with this simple case, just to understand 
qualitatively how ChPT can also account for the description of DCC's. A better 
approximation would be to solve the coupled equations for the o and it fields, 



which yields the solution for a(t) in terms of elliptic functions [15]. Therefore, in 
this simple picture, we take our f(t) of the same form as the lowest order a(t) in 
the LSM, i.e, 
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Here, q is a small parameter, playing the role of <7o in the LSM. Notice that our 
nonequilibrium chiral power counting demands qM 2 = 0(p 2 ) and so on. Thus, for 
defmiteness, we will take q = 0(p 2 /A 2 ), so that the 0(p 4 ) corrections remain under 
control (see below), and M arbitrary. In the end, we will discuss how the results are 
affected by T h q and M. Therefore, we have m 2 (t) = -(qM 2 /2) cosMt(l + 0(q)), 
so that the differential equation (10) becomes to leading order 



d 2 4fc 2 o 

- + -- 2g cosMt 
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where we have Fourier transformed in the spatial coordinates only (k 2 = \k\ 2 ). 
The above equation is nothing but the Mathieu equation, which has several well- 
known interesting properties [26,27]. Among them, it admits unstable solutions 
exponentially growing in time, for certain values of 4k 2 jM 2 . This is the simplest 
version of the parametric resonance mechanism. In particular, the instabilities 
develop in bands in k, centered at k n = nM/2, of width Ak n = 0(q n ). Hence, in the 




FIGURE 1. iMGo(k,t) for T, = M and q = 0.1. Spatial momentum and time are measured 
in units of M and M~ l respectively. The instability band for this case lies roughly between 
QAM <k< 0.6M. 



approximation we are working, only the first band is relevant, i.e, unstable solutions 
only exist for Mjl — Ak± < k < Mjl + Ak\. This is known in the Cosmology 
literature as the narrow resonance approximation [17]. A typical unstable solution 
Go(k,t) has been plotted in Figure 1 for a particular choice of the parameters 
in the first band. The solutions typically oscillate with an exponentially growing 
amplitude inside the unstable region. Therefore, we see that our ChPT approach 
allows for DCC-type configurations. 

Next, we will apply our results for /„-(£) to this particular case. The equal- 
time correlation function G (t) = J d 3 kG (k,t,t) turns out to be UV divergent, as 
expected. After standard manipulations in dimensional regularisation (d — 4 — e) 
one can cast the divergent part for t > as 

t Gf^t) = -^ cos Mt(- + -log -^--] (15) 

v ; 16tt 2 \e 2 M 2 J y ! 

which is an example of the new time-dependent divergences we were talking about 
in previous sections. In fact, we see that it has exactly the same form as f(t)/f(t). 
Furthermore, replacing in (8)- (9), we find that the result is rendered finite and scale 
independent with the same renormalisation of L n derived in [22]. 

The final results for f%(t) are plotted in Figure 2 for different choices of the 
parameters. We clearly observe the damping effect on the amplitude due to the 
unstable solutions at long times. In other words, the DCC's accelerate thermalisa- 
tion. We also observe that this mechanism becomes less efficient for smaller q and 
M. Typically, the unstable corrections to the amplitude of f n (t) are proportional 
to (qM 2 /47rf%) exp(gMt). On the other hand, this effect seems to be rather in- 
sensitive to the initial temperature and thus we expect to catch all the important 
qualitative behaviour concerning the DCC's, regardless of the initial conditions. It 
should be pointed out that the curves have been cut off at the times where the 
one-loop contribution becomes of the same size as the tree level one. From that 
point onwards, the exponentially growing correlator dominates, yielding unphysical 
results. As commented above, we do not expect our simple cosine shape for f(t) to 
be valid for all times, since it is derived neglecting the pion correlator. This final 
time tf roughly defines the applicability range of our results. We expect that this 
range is enough to account for all the plasma time evolution of a realistic RHIC. 
For instance, for M— 1 GeV and q = 0.1, we get tf ~ 35M -1 ~ 8 fm/c. This is 
exactly the same as extrapolating the equilibrium result (11) to predict the critical 
temperature at T = T c ~ 6f%, where all the higher order corrections become of 
the same order. Nonetheless, that formula predicts the right behaviour of f w {T) as 
it approaches the transition. In the same way, our results reproduce the expected 
qualitative behaviour as we extrapolate them up to times t ~ tf. Therefore, /„-(£) 
can be regarded as an alternative observable (it is the residue of an axial-axial cor- 
relator and it can be measured in semileptonic decays) to test the size of DCC-like 
configurations in the late stage of the plasma expansion. 
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FIGURE 2. Results for /^(i) in the parametric resonance regime. The horizontal line is the 
value of /„■ at the initial temperature Tj 



CONCLUSIONS AND OUTLOOK 

We have reviewed recent work on the extension of ChPT to a nonequilibrium situ- 
ation. The NLSM with a time dependent pion decay constant provides a nonequilib- 
rium effective model with a well-defined perturbative expansion and power counting 
near equilibrium. The analogy of this model with curved space-time QFT allows 
to consistently construct higher order lagrangians and implement renormalisation. 
As a first application, we have obtained the renormalised one-loop f n (t). 

We have also shown how this model can be applied to describe DCC-like struc- 
tures in the late stage of the expansion of a hot plasma formed after a RHIC. Work 
in progress includes a more realistic study of the parametric resonances, including 
consistently the pion correlations in f(t) and calculating the correlation length and 
the number of pions. One can also think of including pion masses, extending the 
results to three flavours, using large N methods and cosmological applications as 
other interesting aspects to be investigated in this context. 
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